The distribution of annihilation luminosities in dark matter substructure 
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We calculate the probability distribution function (PDF) of the expected annihilation luminosities 
of dark matter subhalos as a function of subhalo mass and distance from the Galactic center using 
a semi-analytical model of halo evolution. We find that the PDF of luminosities is relatively broad, 
exhibiting a spread of as much as an order of magnitude at fixed subhalo mass and halo-centric 
distance. The luminosity PDF allows for simple construction of mock samples of 7-ray luminous 
subhalos and assessment of the variance in among predicted 7-ray signals from dark matter annihila- 
tion. Other applications include quantifying the variance among the expected luminosities of dwarf 
spheroidal galaxies, assessing the level at which dark matter annihilation can be a contaminant in 
the expected 7-ray signal from other astrophysical sources, as well as estimating the level at which 
nearby subhalos can contribute to the antimatter flux. 

PACS numbers: 95.35.-|-d,98.35.Gi, 98.35.Pr, 98.62.Gq 



I. INTRODUCTION 

Multiple lines of observational evidence have estab- 
lished the existence of a form of non-baryonic dark matter 
binding galaxies. Dark matter is commonly considered 
a yet-to-be discovered elementary particle. A promis- 
ing candidate is a Weakly-Interacting Massive Particle 
(WIMP) that arises in extensions to the standard model 
of particle physics. Examples include the lightest super- 
symmetric particle [ll, 0] and particle excitations in theo- 
ries of Universal Extra Dimensions Q . WIMPs interact 
via the Weak interaction and were in thermal equilibrium 
in the early Universe. If the dark matter is a thermal 
relic WIMP, the WIMP annihilation cross section can be 
constrained by requiring that the present dark matter 
density is flc ~ 0.23 [J]. Implied values for this cross 
section are of order (av) ^ 10~^^ cm^/s, and candidate 
WIMPs can annihilate to many observable states, such 
as 7-rays and high-energy neutrinos. 

The Cold Dark Matter (CDM) model of cosmologi- 
cal structure formation predicts that dark matter is dis- 
tributed in "halos" in a hierarchical fashion. The host 
halo of the Milky Way is expected to have mass of 
M « h-^MQ d g. Within the CDM model such 
halos are also expected to contain numerous smaller dark 
matter "subhalos" which, in turn, contain subhalos of 
their own, perhaps with masses all the way down to the 
cutoff scale in the primordial density fluctuation power 
spectrum 

A well-studied avenue for possible identification of the 
dark matter is to search for unique products of dark 
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matter annihilation at the center of our Galactic halo, 
where densities are highest [TTI-[l8l|. However, astrophys- 
ical backgrounds are also highest toward the Galactic 
center, so an alternative approach is to search for an- 
nihilation within subhal os 11 9143 ill (including potentially 
sub-solar mass halos [sol. I32| - l36j l) . or to infer annihilation 
within subhalos statistically through the angular distri- 
bution of diffuse 7-ray emission |37H47| . More recently, 
the one-point 7-ray flux probability distribution function 
(PDF) was proposed as another way to quantify the ex- 
pected annihilation signal [4^ . The main idea behind 
this approach is that the pixel-to-pixel flux variation from 
7-ray emission from small dark matter subhalos would 
deviate from the Poisson fluctuations that would be ex- 
pected from a smoothly-distributed dark matter halo and 
smoothly-distributed backgrounds. 

A calculation of the contribution of annihilation prod- 
ucts from subhalos to the flux along any sight line con- 
sists of the following two ingredients. The flrst is the 
number of subhalos intercepted along the line of sight, 
assuming subhalos to be small compared to the angular 
resolution of the instrument as is the case with contempo- 
rary detectors. The second is the annihilation luminosity 
of each intercepted subhalo. The latter depends on the 
distribution of dark matter within subhalos, which re- 
flects the underlying process of nonlinear mass assembly 
individual to each subhalo. Although two subhalos may 
have the same mass and be located at the same Galacto- 
centric distance today, they may have different annihila- 
tion luminosities because they may have different forma- 
tion times, mass assembly histories, and different orbits 
in the Milky Way potential. Each of these factors affects 
the internal densities of subhalos. Therefore, we expect 
a distribution of luminosities at each subhalo mass and 
Galacto-centric distance. 

In this paper, we estimate the probability distribution 
function (PDF) of the subhalo luminosity as a function 
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of the subhalo mass and Galacto-centric distance from a 
large ensemble of subhalo populations generated using a 
semi-analytic model of halo and subhalo evolution [48l |. 
Our aim is to use a statistically-large sample of subhalo 
properties to provide a useful tool for computing sub- 
halo annihilation signals in a way that captures some of 
the complexity of nonlinear evolution. We show that the 
luminosity PDF of subhalos is well fit by a log-normal 
distribution, reflecting the underlying distribution in for- 
mation times (or concentrations, pg . Isoj ) and provide 
simple, empirical fits to the luminosity PDFs as a func- 
tion of mass and distance. Applications for the derived 
PDF range from quantifying the variance in the expected 
luminosities of dwarf spheroidal galaxies to generating 
mock synthetic 7-ray sky maps, to understanding the 
level at which dark matter annihilation can be a contam- 
inant in the expected 7-ray signal from other astrophys- 
ical sources j5l|-|54|. and to ascertain the level at which 
a nearby subhalo can contribute to the measured flux of 
antimatter 15511 . 



II. THE IMPORTANCE OF THE SUBHALO 
LUMINOSITY PDF 

The effects of a distribution of luminosities on any cal- 
culation that involves the diffuse emission is only impor- 
tant if most of the diffuse flux originates from a large 
number of sources. If, for example, the diffuse flux is 
due to very few sources with high luminosities, then the 
properties of the diffuse background would be dominated 
by the Poisson statistics of the emission of photons from 
these sources. If, on the other hand, there are many dim 
sources along the line of sight, the intrinsic variation in 
the luminosities of these sources will have an effect on 
the flux PDF. This is due to the fact that the flux PDF 
will deviate from Poisson statistics as it will depend not 
only on the flux-density distribution but also on the mean 
number of sources (see the ^P{D) analysis' discussion in 
the Appendix of Lee et al. ) . 

We can quantify this argument as follows. Suppose the 
number density of subhalos of luminosity L at position £ 
along a particular line of sight is given by dN/dMdV oc 
M~°'n[£{r)], where n[£{r)] is the number density of sub- 
halos at Galacto-centric distance r along line of sight dis- 
tance £. Additionally, assume a mapping between the 
luminosity and the mass of a subhalo, L (x W^. Then 
the contribution to the received flux that is produced by 
subhalos in a given logarithmic mass and line of sight 
interval is given by 



dF 



din Mdlni 



oc M^"+^+i^n[r(£)], (1) 



Taking a « 1.9 and /3 ~ 0.8 (consistent with analytical 
arguments [1^, and numerical simulations [26, 30, 56]), 
the flux per logarithmic interval in mass and line of sight 
distance has a weak dependence on mass Fmi oc Af"'^'^. 



In order for the low mass subhalo contribution to the 
annihilation flux to be roughly the same as high mass 
halos, their abundance must be larger. The mean num- 
ber of subhalos per logarithmic line of sight interval and 
logarithmic mass interval is 



dN 



dhvMd\\\> 



cx M-"+i^3n[r(£)]. 



(2) 



Assuming that the distribution of subhalos traces that of 
dark matter, n[f(£)] cx [f(l -I- f)]~^ [13, Ull, where f = 
r/rg and is the scale radius of the NFW dark matter 
profile. For a w 1.9 and small distances (^ ^ r ^ Ts), 
n[f{l)] cx l/f{l) so that Nmj oc M-°-^£'^. For large 
distances {£ fa r :s> r^), n[r{£)] cx l/f^{£) and Nm,£ oc 
]\,f-o.9^ roughly independent of £ for all masses. 

The previous two paragraphs show that we should ex- 
pect the majority of the diffuse annihilation fiux to be 
due to the presence of numerous low-mass subhalos. In 
the left panel of Fig. [T] we show the quantity Fm/. in units 
of cni~^s~^sr~^, as a function of subhalo mass and posi- 
tion along the line of sight, Eq. ([T|), at an angle ip = 90° 
with respect to the Galactic center. We assume a Milky 
Way halo with a radius Rmw = 250 kpc, and a scale 
radius « 20 kpc. Note that at any fixed mass, most 
of the fiux comes from a region at £ « 10 kpc. This 
is because the number of objects declines rapidly with 
distance at Galacto-centric distances signflcantly larger 
than a scale radius. Moreover, low-mass subhalos con- 
tribute marginally more flux than their high-mass coun- 
terparts. In the right panel of Fig. [T]we show the mean 
number per logarithmic mass and line-of-sight interval, 
Nm,i- The two panels of Fig. [T] illustrate a basic conclu- 
sion that the mean number of objects along a line of sight 
increases with decreasing subhalo mass (right panel) , and 
as all intervals of subhalo mass have comparable flux con- 
tributions (left panel) , the signal is set by relatively low- 
mass subhalos that are close to the observer. This is in 
qualitative agreement with the result of Lee et al. [i^ 
who found that substructures give rise to photon counts 
that deviate from a Poisson distribution. 

Our simple demonstration neglects the presence of a 
baryonic disk in the Milky Way halo, and its effect on 
the subhalo population in the inner regions of the halo. 
Recent studies flnd that the inner 30 kpc of a Milky Way- 
sized halo may be deflcient in subhalos due to interactions 
with the disk [5^ . This should have an effect on the flux 
contributions we derived, but the approximate mass and 
distance dependence that we describe should be main- 
tained. Nevertheless, it is important to keep in mind 
that a depleted subhalo population in the inner regions 
of the Galactic disk may have an effect on the fiux from 
substructure. The maximal net suppression in annihila- 
tion flux near the disk due to this suppression is expected 
to be a factor of a few ISOll. 
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FIG. 1: Left: The mean contribution to the annihilation^ flux from all subhalos in a logarithmic interval of subhalo mass 
(horizontal axis) and line-of-sight distance (vertical axis), FM,e defined in Eq. ([1} in units of cm~^s~^sr~^. This panel shows 
the mean contributions along a line of sight 90° away from the Galactic center. The color at each position designates the value 
of Fm,i, where lighter colors indicate larger values as shown in the color bar. The contribution along the line of sight peaks 
at ~ 10 kpc at all masses. Right: The mean number of subhalos per logarithmic interval of subhalo mass and logarithmic 
line-of-sight interval. Lighter colors designate larger numbers of subhalos according to the color bar. Smaller objects dominate 
the number density at all radii [see Eq. ((2])]. 



III. THE SUBHALO ANNIHILATION 
LUMINOSITY PDF 

In order to derive the subhalo luminosity PDF, we 
model the accretion and dynamical evolution of subha- 
los in the Milky Way potential well using an approxi- 
mate semi-analytic technique [i^ . This approximate ap- 
proach greatly reduces the computational cost of model- 
ing substructure by treating subhalo density profiles as 
continuously-evolving functions that can be described by 
a small number of parameters, rather than a collection of 
a very large number of individual particles. This enables 
calculations of the properties of subhalo populations in 
a large set of Milky Way-sized halos, so that object-to- 
object variance can be estimated. This is not yet feasible 
in simulations directly. This procedure also provides a 
physically-motivated method to extrapolate the results 
of numerical simulations to masses below their resolution 
limits, masses which are not negligible from the annihi- 
lation perspective. Of course, the cost of this method 
is that it is approximate and the approximations imple- 
mented cannot be validated outside the range of scales 
that are resolved by A^-body simulations. This method 
is described in detail in , which also shows a number 
of non-trivial comparisons with A^-body simulations that 
validate this treatment of halo substructure. 

We use an ensemble of 200 realizations of subhalo pop- 
ulations within a dark matter halo of mass A/mw = 
1.26 X lO"'^^/i~^M0 in a flat cosmological model with 
= 0.228, ^bh^ = 0.0227, h = 0.71, and erg = 0.81, 
favored by the five-year Wilkinson Microwave Anisotropy 



Probe results f3J. Each realization represents a possible 
subhalo population within a Milky Way-like halo. The 
populations differ because only some statistical proper- 
ties of the initial density field in the local neighborhood 
are known, not the precise initial conditions for collapse. 
We use this large number of distinct populations to quan- 
tify the predicted variation from one Milky Way-sized 
halo to another. The result of the calculation is a list of 
all subhalos, complete with their structural parameters 
such as bound mass, scale radius, and tidal radius. These 
quantities all evolve with time, and we study them at the 
present epoch (z = 0). 

Figure [2] shows the cumulative velocity function of 
Galactic subhalos. The power-law behavior is similar to 
what is found in numerical simulations, with the number 
of subhalos increasing as A^(> T4iax) « V^^^. Our model 
also reproduces the abundance of subhalos in cosmolog- 
ical simulations reasonably well (see also Ref. for a 
discussion) . For example, the average number of subhalos 
with Umax > 4 km/s within the inner 200 kpc of the Milky 
Way in our models is 2481 with a 68 percentile range 
of [1964-3007]. This is consistent with 2469 subhalos of 
I'max > 4 km/s within the same radius in the Via Lactea 
II simulation of a Milky Way-sized halo [g^l- Moreover, 
the dispersion in the number of subhalos is consistent 
with the empirical finding of a Poisson scatter added in 
quadrature to an intrinsic scatter of about 20% derived 
from recent numerical simulations [6l|. We note that 
the validation exercises in Ref. [i^ show that this model 
may over-estimate halo-to-halo variance mildly (see their 
Fig. 7), though it is difficult to assess through a compari- 
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son with simulations with greater precision at this point. 

Note that in general, the slope and normalization of the 
velocity function depend on the concentration of the host 
dark matter halo [48]. As shown in Fig. [21 the normal- 
ization of the velocity function of subhalos is reduced. 
This is to be expected as host dark matter halos with 
high concentrations are formed earlier than low mass ha- 
los. As a result, there is more available time for evolu- 
tion of the subhalo population leading to a decrease in 
the total number of subhalos. A subsidiary effect that 
also contributes to the decrease of the subhalo popula- 
tion of high-concentration halos is the efficient tidal dis- 
ruption due to the higher central densities. On the other 
hand, subhalos in host halos with low concentrations will, 
on average, have spent less time in the host and have a 
higher rate of survival due to the decreased tidal forces. 

We model the final dark matter distribution in a sub- 
halo of mass M as a Navarro- Frenk- White (NEW) profile 



p(r) = 



M 



1 



1 



A-Krl f{ft) r{l + fy 



(3) 



In Eq. ([3]), f = r/vg, where is the evolved scale ra- 
dius. f{x) = ln(l -I- a;) — x/(l -I- x), and ft — rt/rs is the 
ratio of a tidal truncation radius r^, to the scale radius. 
We treat tidal truncation by assuming an abrupt limit 
to the extent of the subhalo density profile for simplicity. 
Isolated high-resolution simulations of subhalo evolution 
support a very sharp truncation (e.g., Ref. [6^, [63|) and, 
furthermore, the annihilation luminosity emanating from 
the transition region should be very small in comparison 
to the total annihilation luminosity. Our simplifying as- 
sumption of a truncated NEW profile is not crucial to our 
analysis. Deviations in the assumed power law in the in- 
ner regions of the profile do not result in large changes 
in the total luminosity of the halo. Profiles with differ- 
ent inner slopes also require different normalizations to 
maintain the same bound mass against the tidal forces 
of the host. Moreover, much of the total luminosity of 
a nearly-NEW halo arises from the region near the scale 
radius, r^. As such, changes in the inner slope do not 
lead to large changes in the annihilation luminosity in 
most circumstances (see Refs. [23l. [27l[29l[6^ |65| ) . 

The dark matter luminosity of a subhalo is then ob- 
tained from 



L = 47r 



Ml 



2 2 J 

p r ar 



3.32 X lO^^ph {(jv)-2&N: 



7,30 



kpc/ 



p{7 



GeV/cm 



r^df. 



(4) 



Here, ((Tf)_26 is the annihilation cross section in units of 

is the mass of the dark matter 
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FIG. 2; The cumulative velocity function of subhalos. The 
solid black line shows the average velocity function over all 200 
realizations of the formation of a Milky Way-size halo. The 
short-dashed red line shows the velocity function derived for 
Galactic halos which have a concentration c > 13.4, while the 
long-dashed blue line shows the velocity function for Galactic 
halos with c < 6.7. The dotted line shows the behavior of a 
A'^(> Vmax) oc K~ax power law, which describes the subhalo 
velocity function in cosmological simulations. 



units of 30. This fiducial choice of parameters is represen- 
tative of optimistic scenarios in the Minimal Supersym- 
metric Model. We emphasize that we defined luminosity 
as number of photons per unit time, and not energy per 
time. Similarly, when we discuss flux, we implicitly mean 
photon flux, and not the energy flux. Note also that as 
the goal of this work is to quantify the spread in lumi- 
nosities in the subhalo population, the choice of particle 
physics parameters in Eq. |4] is intended only to provide a 
useful representation of the magnitude of the luminosity. 
However, the PDEs of luminosities derived below are af- 
fected by the choice of particle physics parameters only 
in their normalization. 



IV. RESULTS 

We compute subhalo luminosities as a function of ra- 
dial position and mass as follows. We flrst determine the 
minimum and maximum of the radial distribution (typ- 
ically Tmin ^ kpc, Tniax ^ 250kpc) of subhalos as well as 
the minimum and maximum subhalo mass at the present 
epoch (typically M^r^m IO^Mq, M,„ax 10"Mq) . 
We then divide the radial distribution of subhalos into 
Nr bins. Each radial bin is then subdivided into Nm 
mass bins of equal logarithmic size in mass. We then flt 
the distribution of luminosities for different values of Nr 
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FIG. 3: Example probability distribution functions of the 7-ray luminosity for three sample mass and radial bins. The solid 
curve depicts a histogram of the distribution of luminosities in each bin. The dashed curve shows the fitting function of Eq. [5] 



and Nm until the maximum deviation of the fit is less 
than 10% of the true value. We find that N„i ^ Nr ^ 50 
with at least 200 subhalos per bin has errors of at most 
10% of the true value. 

The luminosity of a subhalo is a function of its mass 
as well as its position within the host halo. We find that 
the subhalo luminosity PDF is well fit by a log-normal 
distribution, as 



P[\nLM,r] 

where, 



1 



1 
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and, 



O'ALr = 61 + ^2 In 



M 



(io^Mq 



63 In 



50kpc 



(5) 



(6) 



(7) 



Here, we implicitly assume that the luminosity LM,r is 
expressed in units of photons per second. 

Assuming the fiducial particle physics parameters as 
shown in Eq. 21 the best fit parameters for the whole 
population of subhalos in all 200 realizations are ai — 
77.5, 02 = 0.87, 03 = -0.22, 61 = 0.75, 62 = -0.0026, 
and &3 = 0.0061 (see "All" in Table U). This resuh can 
be scaled to any assumed particle physics parameters, by 
simply adding the term 



In 



9 X 10~29cm3s-iGeV-^ 



(8) 



to the parameter oi of Eq. [HI 

The fitting function [Eq. [5] with Eq. © & Eq. ©] is 
good to within ^ a few % for P[\nLM,r] as a function of 
In Lm^t over the range of masses and radii we have exam- 
ined with sufficient statistics, M sa [10^ — 10^°] Mq, and 





ai 


a2 


ffls 


bi 


62 


63 


All 


77.4 


0.87 


-0.22 


0.75 


-0.0026 


0.0061 


Co 


77.4 


0.87 


-0.23 


0.74 


-0.0030 


0.011 


c+ 


77.5 


0.87 


-0.26 


0.76 


-0.0021 


0.0077 


C- 


77.3 


0.87 


-0.18 


0.75 


-0.0013 


0.0043 



TABLE I: Fitting parameters for the mean and width of the 
7-ray annihilation flux distribution function (see text). 



r sa [5 — 250] kpc. The mean and variance of the distribu- 
tion are functions of mass and radius, refiecting the fact 
that the PDF of luminosity is set by the interplay be- 
tween the mass function of accreted objects, the redshift 
of accretion, and the orbital evolution of the individual 
subhalos constituting the population. 

The mean luminosity of isolated, field halos that have 
not experienced strong interactions within the potential 
of a more massive halo, scales as L cx p^r^ oc Mc^ / p{c), 
where M is the halo mass, corresponding to a virial ra- 
dius radius R cx M^/^, c = R/vg is the concentration and 
/(c ) (X c°'* near c « 30 as is relevant for small subhalos 
(27|. If we assume a weak dependence of concentration 
on mass c(M) cx M-° i [H, l50|, l66l468|. the luminosity 
of a halo scales roughly as iy oc M°-^. However, sub- 
halo populations deviate from this scaling somewhat for 
several reasons. Subhalos merge into the host halo at a 
variety of times, so they sample the c{M) relation at a 
variety of redshifts and subhalos of fixed mass at the time 
of merger exhibit a variety of bound masses at z = as a 
result of their distinct orbital evolution histories. There- 
fore, subhalo luminosities at fixed mass are infiuenced by 
the mass and redshift dependence of concentrations and 
subhalo orbital properties. The effects of these physical 
changes in the structure and description of subhalos com- 
pared to isolated halos are reflected in the best fit param- 
eters of Eq. [6] In particular, we find that the luminosity 
of evolved subhalos scales approximately as L oc M°-*^. 

In Figure |3] we show the distribution of luminosities 
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FIG. 4: Probability distribution functions of the 7-ray luminosity of a 10 subhalo at three different radial distances from 
the Galactic center (Eq.[5]and parameters from Table H]). The solid curve depicts a histogram of the distribution of luminosities 
in host halos with a concentration in the 68 percentile about the mean of the sample of 200 realizations (Co). The short-dashed 
curve shows the subhalo luminosity PDF for host halos with a concentration in the upper 16 percent (C+), while the long-dashed 
curve shows the subhalo luminosity PDF for host halos with a concentration in the lower 16 percent (C_). 



in three examples of radial and mass bins. The figure 
also shows the fit to the luminosity distribution given 
in Eq. ([5]). The choice of mass and radius displayed in 
Fig. [3] are only meant to demonstrate schematically the 
agreement between the fitting functions and the numeri- 
cal results. Notice that the spread in luminosities at fixed 
mass and position can be roughly an order of magnitude, 
and that the peak of the luminosity PDF depends on 
both, mass as well as radius. 

As the structural properties and abundances of subha- 
los are influenced by the merger history and properties 
of the host halo, we expect that the luminosity PDF will 
be affected by the distribution of dark matter within the 
host dark matter halo. In order to explore the relation- 
ship between the host halo properties and the subhalo 
luminosity PDF we divided the subhalo populations into 
three groups based upon the concentrations of their host 
halos. We grouped subhalos with host halo concentra- 
tions within the 68% range of the mean concentration, 
c ~ 9.7. This group, "Co", consists of subhalos in pri- 
maries with concentrations in the range 6.7 < c < 13.4. 
We then collected all subhalos in hosts in the upper and 
lower 16% ranges, with group C_ consisting of host con- 
centrations in the range c < 6.7 and group C+ consisting 
of primaries with c > 13.4. The fitting parameters of 
Eq. ([U for these three concentration bins are also given 
in Table m 

Fig. U shows the effects of the host halo concentration 
on the luminosity PDF of subhalos. At inner radii (left 
panel of Fig. S]), the effects are more pronounced. This is 
the region where earliest merging subhalos reside. Host 
halos with high concentrations are formed early, and con- 
tain subhalos that on average have formed earlier, and 
so are also more concentrated [13]. As the annihilation 
signal is sensitive to the concentration of subhalos p7| . 
the luminosity of subhalos in a high concentration host 



is slightly higher (C+ curve in Fig. |4]) than subhalos re- 
siding in host halos with lower concentration (C_ curve 
in Fig. [5]). At large radii, r > (right panel of Fig. 2]) 
the difference between host halos of different concentra- 
tions diminishes, reflecting the fact that the outer regions 
of halos typically contain recently-merged substructure 
(the host concentration influences subhalo dynamics lit- 
tle beyond the halo scale radius). Fig. 2] and Table |T] 
demonstrate that the maximum shift in typical subhalo 
luminosities is in the inner regions of host halos, and is 
relatively small, ~ 20%. 



V. APPLICATIONS 

The flux PDFs that we provide have numerous po- 
tential applications. As an example, in this section we 
compute the probability distribution for measuring a flux 
between F and F + dF from any individual subhalo along 
a line of sight at an angle tp, from the Galactic center. 
We compute this "single halo flux PDF" as 



Piim 



oc 



e 



4 dN[r{i,^)] 



J M„i„ dMdV 
P[LFj.\M,r{e,i:)]dM di 



(9) 



Here, £ is the line-of-sight distance, dN[r{£,i/j)]/dMdV 
is the subhalo mass function, and Lpj — AttI'^F ensures 
a proper flux measurement for a subhalo of luminosity 
L at a distance I. The quantity £max is the maximum 
line-of-sight distance we consider, given by 



= dr. 



cos V' + \l {Rc/dey - sin^ ip 



(10) 



where dQ — 8 kpc is the distance of the Sun to the Galac- 
tic center, and i?G = 250 kpc is the approximate radius 
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of the Galactic halo. As in Sec. |TT1 we assume a mass 
function of the form dN/dMdV cx M-"/f(l + f^With 
a = 1.9, as predicted by both N-body simulations [69|, as 
well as by the semi-analytic model of subhalo populations 
we use for this study. Both simulations and our analytic 
method show little evidence that a varies significantly as 
a function of radius [sOjIH,!!^, so any variations are sub- 
tle (though they may depend upon global host halo prop- 
erties) . Consequently, it is convenient and informative to 
couple standard halo mass functions with our luminosity 
PDFs to estimate relevant observable quantities. How- 
ever, we emphasize that this particular choice of mass 
function is not unique and, in principle, the flux PDF 
can be derived using Eq. ([5]) with a variant subhalo mass 
function. 

In Figure [5] we show the normalized single halo flux 
probability distribution function for two different lines 
of sight generated from the contributions of subhalos 
in the range Mmin = [10"' - lO^°]Af0 (thick lines) and 
Mmin = [10^ - 10^°]Mq (thin lines). This choice of sub- 
halo mass does not affect the overall shape of the mass 
function. A change in the value of Mmax does not affect 
the result as the mass function power law is a steep de- 
creasing function of mass (see also the left panel of Fig. [l] 
which shows how the average flux due to high-mass halos 
is lower than that due to low-mass halos). The high- flux 
power-law shape of the flux PDF can be understood in 
the following way. As the mass function is proportional 
to dN/dMdV cx M~^-^, and the luminosity of a subhalo 
scales with mass as L cx A/°-*^ (see Table I), then the in- 
tegrand is proportional to L^^ '^^. However, F cx L/£'^, so 
the flux PDF is a power law with a shape given roughly 
by P[F] cx rpj^jg jg apparent for both lines of 

sight in the high-flux regime. 

Changes in the value of M^nin affect the low-flux be- 
havior of the flux PDF. It is easier to understand the 
low-flux cutoff if we assume a delta function lumin osity 
PDF instead of the log-normal distribution, as in [43| . 
For a line of sight at some angle tp, there is a maxi- 
mum distance ^max which is a function of '0 and the ra- 
dius of the Milky Way halo. If we assume that L cx Af , 
and that the subhalo radial distribution is the same for 
all masses, then the minimum flux would be given by 
i^min c>c Afmin/^maxi i-^-: the Smaller the minimum mass, 
the broader the PiiF). For illustrative purposes we show 
the Pi{F) flux PDF derived under the assumption of a 
delta function luminosity PDF as in [43] in Fig. [5j As- 
suming the log- normal luminosity PDF (instead of a delta 
function), results in a tail in the low- flux region of the 
Pi{F), and thus a broader Pi{F) PDF. 

The slope of the flux PDF for low fluxes does depend 
on the angle between the line of sight and the Galactic 
center (see the blue solid, and red dashed lines in Fig.[5|). 
At small angles, there is flux probability excess in the 
low flux regime relative to the high flux region. This 
can be explained using Fig. [T] At large angles from the 
Galactic center, the total number of halos intercepted 
along a line of sight is smaller than the total number of 
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FIG. 5: The flux probability distribution function derived 
using Eq. [5] and the Co parameters of Table [D The solid 
blue lines is the flux PDF along a line of sight at 20° relative 
to the Galactic center, while the dashed red lines are along 
180°. Thick lines correspond to a minimum subhalo mass of 
lO^M©, while thin lines correspond to a minimum subhalo 
mass of 10^ Mq. The black dotted line depicts the single halo 
flux PDF derived under the assumption of a delta function 
luminosity PDF. 

halos intercepted along a line of sight that passes near the 
Galactic center. In addition, the fraction of these halos 
that are close to the Sun is smaller for large ip. Therefore, 
any effects due to the spread of luminosities will be more 
pronounced where the fraction of contributing sources is 
large, and also nearby (due to the l/£^ term). The spread 
of luminosities thus introduces a spread in flux, giving 
rise to substantial a change in the power-law behavior of 
the flux PDF. 

It is important to draw the distinction between the 
single halo flux PDF (Eq. [3]) and the probability distri- 
bution function of measuring a total flux P{F), from the 
contribution of numerous halos along the line of sight. 
The total flux PDF can be computed in a straightfor- 
ward manner from the basic quantity Pi (F) [70-74] , and 
more recently [i^). A thorough investigation of the total 
flux PDF using the single halo PDF Pi (F) , is presented 
in 75] , where the authors evaluate the ability of FGST to 
discover dark matter via 7-rays from Galactic substruc- 
ture. 

Another possible application of the luminosity PDF 
is in the intrinsic spread of the expected 7-ray luminos- 
ity signal from dwarf spheroidal galaxies in the Milky 
Way halo. Dwarf spheroidals are very low surface bright- 
ness systems dynamically bound to the dark matter halo 
of the Milky Way. Their very high mass-to-light ratios 
make them ideal for 7-ray studies as their dark mat- 
ter distributions can be well constrained using the ve- 
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locity dispersions of their stellar populations. Numerous 
studies addressed the possibility of detecting dark mat- 
ter annihilation in dwarf spheroidals using either FGST 
or ground-based Cerenkov telescopes, such as VERITAS 
and H.E.S.S. ^19ri23., .31] . Dynamical studies of the ve- 
locity dispersion of stars in these systems constrain the 
dark matter mass (and profile) [2^, [3l[- Suppose that 
a group of dwarf spheroidal galaxies have masses es- 
timated via dynamical measurements. The luminosity 
PDF presented in Sec. |IV] can be used to assess the ex- 
pected variance of a possible detection of 7-rays from 
either Cerenkov telescopes or the FGST. The estimated 
variance can be used in concert with the measurement of 
the distribution of dark matter from dynamical studies 
in the interpretation of any signals or limits from these 
systems. 

Moreover, the use of the luminosity PDF can be impor- 
tant in studies aimed at disentangling the different source 
contributions to the diffuse gamma-ray background, such 
as blazars, starburst galaxies, pulsars, supernova rem- 
nants, and potentially cataclysmic binary systems, and 
dark matter [38-42, 4^ \M, MM,^- A useful tool in 
these studies is the use of the angular correlation function 
of flux fluctuations, which is estimated from the convo- 
lution of the emissivity as a function of distance and the 
spatial correlation function of the sources. Knowledge of 
the expected distribution in luminosities of contributing 
dark matter substructure can be used to remove the dark 
matter contamination to the astrophysical background, 
and thus potentially open the window for the detection 
of 7-rays from yet undiscovered sources [H, [t^, [z1| • 

Finally, the annihilation luminosity PDF can be used 
in studies aimed at estimating the likelihood that nearby 
dark matter subhalos contribute significantly to the mea- 
sured antimatter flux [t^Is^I- A dark matter explanation 
for the anomalous excess in antimatter flux must rely on 
either non-standard extensions to the standard model of 
particle phy sics (prominent annihilation to charged lep- 
tons j8ll - [84| . or some new, long-range force 85], or po- 
tentially the presence of a dark matter halo in the near 
solar-system neighborhood ^]. The luminosity PDF 
presented in this paper can be used to assess the min- 
imum spread of the expected flux from a given subhalo 
and the likelihoods of particular flux measurements from 
rare, nearby objects. Note however, that propagation ef- 
fects as well as the degeneracy between luminosity and 
the square of the distance to the subhalo, and velocity- 
dependent annihilation [gl, [s^ [s^l can all increase the 
spread considerably. 

VI. CONCLUSIONS 

We have presented an estimate of the luminosity prob- 
ability distribution function of dark matter subhalos 



within a Milky Way-like parent halo. An empirical flt 
to our numerical calculations suggest that the luminos- 
ity PDF can be described well by a log-normal distri- 
bution with subhalo mass- and position-dependent mean 
and variance. This log-normal probability distribution 
has a width that is determined by the wide distribution 
of formation times and concentrations for both the host 
halo and subhalos. 

The derived luminosity PDF can be used as an ingre- 
dient in a number of interesting calculations regarding 
predictions for observable dark matter annihilation prod- 
ucts. This distillation of a complex set of halo properties 
should be particularly useful with the impending data 
from the Fermi Gamma-ray Space Telescope and contin- 
ued advances of ground-based Air Cerenkov Telescopes 
as well as neutrino telescopes and antimatter detection 
instruments. The tool we provide can be used to address 
such observations, including variance in the predicted sig- 
nal, in a relatively simple manner. 

As a straightforward example of the application of the 
PDF, we have estimated the distribution of observed 7- 
ray fluxes along lines of sight as a function of the angu- 
lar separation between the line-of-sight and the Galactic 
center. This may be used as a signature to diagnose unre- 
solved annihilation in a population of Galactic subhalos 
(see [75I, where the authors utilize this distribution to 
estimate the robustness of estimating dark matter prop- 
erties from the diffuse flux measurements of FGST) . Ad- 
ditional applications of the annihilation luminosity PDF 
include estimates in the spread of the angular power spec- 
trum of flux fluctuations as a probe of unresolved sub- 
structure in the Milky Way, as well as the antimatter flux 
distribution from nearby subhalos. These applications 
make the luminosity PDF a useful tool in the analysis 
of forthcoming data in the ongoing effort to identify the 
dark matter. 
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